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ABSTRACT 

A synoptic analysis of the Ionian Sea in June 1980 showed 
eaemmal patterns Of various scale sizes, and in particular, 
a warm core eddy which was comparable in size (ca. 30 km) 
and location with one found in MILOC-68. Spatial correlation 
functions of temperature were anisotropic in the southern 
Bart, with zero-crossings of 30 km in the EW direction and 
40 to 80 km in the NS direction, commensurate with the (first 
mode) baroclinic Rossby radius of deformation of 23 km. There 
was a shallow main sound channel (axis at ca. 100m) and some 
secondary sound channels caused by temperature inversions. 
The spatial patterns of these sound channels were associated 
with synoptic scale and mesoscale oceanic phenomena, which 
in turn were influenced by bathymetry. The temporal patterns 
were influenced by atmospheric forcing. It may be possible 
to use features seen in satellite IR imagery to deduce zones 
with shoaler thermocline depth. With data sampling appropri- 
ate to the scales of atmospheric and oceanic variability, a 
valuable synoptic scale analysis can be derived from XBTs 
collected under operational conditions. Recommendations for 
a regional approach to military oceanography and XBT sampl- 


ing procedures are made. 





RESUME 


L'analyse synoptique de la mer Ionienne en Juin 1980 
révele des structures thermiques variées, et en particulier, 
un tourbillon chaud comparable en taille et position a celui 
trouve au cours de MILOC-68. Les fonctions de corrélation 
Spatiale en température sont anisotropiques dans la partie 
sud, avec des échelles respectives de 30 et 40 a 80 km dans 
les directions Est/Ouest et Nord/Sud, comparables au rayon 
baroclinique de déformation (ler mode) de 23 km. Le chenal 
sonore est situé a environ 100m de profondeur, et des inver- 
sions thermigques créent des chenaux sonores secondaires. La 
distribution spatiale de ces chenaux sonores est liée aux 
phenoméenes océaniques a échelle synoptique et moyenne, eux- 
mémes influencés par la bathyméetrie. Les variations tempo- 
relles sont liées aux effets atmospheriques. I1 semble pos- 
Sible d'utiliser l'imagerie satellite infra-rouge pour 
determiner la profondeur de la thermocline. I1 est conclus 
qu'une analyse synoptique valable peut é6tre effectuée avec 
des sondages bathythermigques préleves en ambiance opération- 
nelle a condition que ce prélevement tienne compte de la 
variabilité océanigque et atmosphérigque. Des recommendations 
sont émises, concernant une approche régicnale de 1l'océanog- 


raphie militaire et les procédures de collecte bathythermigque. 
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I. INTRODUCTION 


A. FOREWORD 

The history of naval warfare has often demonstrated the 
impact of the oceanic environment on combat capabilities. 
Until recently, little was understood of the behavior of the 
ocean. For centuries, the traditional answer had been to 
develop some kind of innate feeling for the environment, as 
successfully demonstrated by mariners or fishermen. In more 
recent years, with the growing sophistication of naval 
ships, there has been a temptation to ignore the environment 
and consider its effects as harmless to a powerful technology, 
or at least secondary. But reality seems to go precisely the 
other way, and the solution, if any, lies through an increased 
knowledge of atmospheric and oceanic processes, which is 
still a challenge to the scientific community. As partial 
answers start to appear, a demand arises to strengthen the 
links between ocean science researchers and the users of the 
sea. Naval officers need to be informed and convinced that 
Seeano!Ogy can provide useful, practical help in a variety of 
activities, ranging from antisubmarine warfare (ASW) to search 
and rescue. 

This informative role will be one of the author's respon- 
Sibilities in his next assignment with the environmental / 


Operational group of the French Navy. With experience in 


ins 





ASW and a newly acquired knowledge in oceanography, he felt 
the need to test and apply it to a practical environmental 
problem in a way that would include a scientific approach to 
data analysis as well as a goal of practical, operational 


eueDuG products. 


pee OBJECTIVES 
This study was, hence, designed as an attempt to: 

- acquire a keen knowledge and understanding of a 
particular oceanic region, especially its dynamical 
processes. 

°" analyse a set of operational data, i.e., data collected 
by naval ships in the course of an exercise in that 
area. 

* yleld operational products in the field of acoustic 
propagation. 

* compare the results of the analyses with other studies. 

- develop and tailor analytical tools, mainly computer 
software, for future use. 
mincdinly, cO evaluate the validity of this kind of 
computer-assisted, but basically manual approach as a 
Support to naval operations. 

Fortunately, a set of data was found in an area of some 
interest to the French Navy, a fact that increased to some 


extent the benefits of this study. 





C. SUMMARY 


zone: 33 to 39N, 14 to 19E, the Western Ionian Sea in 


the Central Mediterranean Basin. 


Period: 8 to 20 June 1980 


Data: 


167 expendable bathythermograms (XBT) taken during 

a US Navy exercise (SHAREM 38). 

14 infrared AVHRR satellite pictures (from the Centre 
de Meteorologie Spatiale, Lannion, France) 

daily wind field analyses (from the Fleet Numerical 
Oceanography Center, Monterey and Meteorologie 
Nationale, Paris, France) 

3 complete analyese (EOTS) of thermal structure by 


the Fleet Numerical Oceanography Center. 


Computational tools: 


the IBM 3033 for the main computations (Computer 
Center, Naval Postgraduate School) 
the Tektronix 4081 for most graphic applications 


(Dept. of Mechanical Engineering) 


OU BoOuUc : 


Thermal analyses from the surface to the 400m level 
Resolution of small scale anomalies 
Spatial and temporal variability 


Acoustic variables 
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be CONTENTS 

Chapter II of this study reviews the climatology of the 
area, with emphasis on the exchange of water masses between 
the Eastern and Western Mediterranean Sea. Chapter III 
describes in some detail the way data were checked for 
quality. Selection of subareas for detailed analysis of 
thermal features is also described. It defines the statisti- 
cal indices (cross-correlation) which are used later, and 
then addresses the processing of "other" data, e.g., 
satellite pictures and wind field. The description of the 
analysed fields commences with Chapter IV. The atmospheric 
conditions and their possible effect on the oceanic thermal 
structure are explained. Chapter V evaluates the different 
processes which might have a contaminating effect on the 
synopticity of the XBT data set, i.e., advection and non- 
synoptic sampling, and concludes with a reasonable assumption 
of synopticity. Chapters VI and VII describe in much detail 
the thermal structure found in the XBT data set, and gives 
some statistical (cross-correlation) parameters. The thermal 
fields are next analysed to yield acoustic variables, as 
discussed in Chapter VIII. 

Two comparisons are made with other analyses: Chapter IX 
evaluates the results, relative to the Fleet Numerical Ocean- 
ography Center (FNOC) analysis, and Chapter X evaluates the 
results in the light of a previous, extensive study made in 


1968 in the same area (MILOC-658). 





Chapter XI summarizes the oceanographic picture, and, 
in conclusion, Chapter XII examines diverse qualitative 
problems related to the use and conduct of oceanographic 


analyses for and by forces at sea. 
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il. "GENERAL FEATURES 


The Western Ionian Sea, a part of the Central Mediter- 
ranean Sea, is bounded to the north by the Italian peninsula, 
to the south by the African continent, to the west by the 
Seeait OL Sicily while its eastern part is open to the 
Levantine Basin (Fig. 1). Its most prominent characteristics 
are described below in terms of bathymetry, meteorology and 


oceanography. 


meee BATHYMETRY 

The Western Ionian Sea, the deepest basin in the 
Mediterranean Sea, has an average depth of 3300m. The main 
boundary is the Maltese Escarpment, a steep barrier which 
rises to the shelf, at an average depth of 200m (Fig. 2). 

The only opening to the western basin lies between Malta 
and the Medina Bank, with a sill depth of 430m. Hence, no 
intermediate water can be exchanged with the Western Mediter- 
ranean basin below this sill. There is no other deep opening 
either north of Malta or in the Gulf of Syrta. East of the 
Escarpment lies a relatively flat area. However, there is 
One notable topographic feature near 35°N, 17° to 18°E: a 
Significant ridge, which rises 1000m over a 10 km horizontal 


distance, and which is about 90 km long. 


Zale 





B. METEOROLOGY 

The atmospheric regime of the Mediterranean Sea is 
characterized by an alternation between the summer and 
winter seasons, the spring and fall seasons being brief 
transitional periods. The winter season commences abruptly 
around October, with the invasion of cold air into the 
Western Mediterranean. Typical storm tracks originate in 
the Gulfs of Lion or Genoa and storms move southeast across 
Italy and Greece. Starting in February, the dominant storm 
tracks originate in the Atlas Mountains and move northeast- 
ward over the area [Brody and Nestor, 1980]. 

The summer season starts in May and lasts until September. 
As a result of a very stable situation where the Azores High 
extends well over Central Europe, this season is hot and 
relatively dry. 

imewune, the period in which this thesis 1s particularly 
concerned, winds are calm about 25% of the time and average 
4 m/s from SW to NW 50% of the time [Naval Weather Service 
Command, 1970]. Sea breeze regimes are common in coastal 
areas. Cloudiness is low (less than 4/8) 803 of the time, 
which accounts for intense solar radiation over the sea, in 
conjunction with low wind stress, and offers exceptionally 


good conditions for passive satellite remote sensing. 
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C. OCEANOGRAPHY 
1. General circulation in the Eastern BaSin 

Because of the hot, dry climate, the Mediterranean 
Sea acts aS a concentration basin for the Atlantic surface 
water. Evaporation in the extreme eastern part greatly 
exceeds both precipitation and river runoff. The water 
budget and the general circulation features of the Mediter- 
ranean Sea are considered to be well established [Lacombe 
and Tchernia, 1960; Morel, 1971]. 

The slope of the sea surface, due to evaporative 
water loss at its eastern end, is the driving force for water 
renewal, with an estimated cycle of 100 years. Atlantic 
Water enters as a near-Surface current through the Strait of 
Gibraltar and flows eastward (Fig. 3). In the summertime, 
the absorption of solar radiation exceeds, by far, the 
radiative and evaporative heat losses. Hence, this surficial 
Atlantic Water, about 100m thick, warms progressively as it 
moves along the Algerian coast. The salinity 1S approximately 
37.5%.. A thermocline is formed around 20 to 40m depth. As 
this tongue of water moves eastward through the Strait of 
Sicily, it splits into two branches. In the Ionian Sea, one 
branch continues eastward, while the other flows southward 
and appears to be modified by the bottom topography [Morel, 
fom Grancini, 1972; Briscoe et al., 1974]. 

Temperature and salinity increase in the surface layer as 


it reaches the Eastern Basin. Below the thermocline the 
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Atlantic Water maintains its relatively low temperature and 
Salinity characteristics, which can be traced to the eastern 
part of the Sea in the Levantine Basin. Hence, because of 
the very strong thermocline, the direct influence of atmos- 
pheric forcing and air-ocean interaction is concentrated in 
a shallow layer. 

In fall and winter the thermal budget in the Eastern 
Basin becomes negative, i.e., the heat flux is from the ocean 
to the atmosphere. Sea-surface temperatures may exceed air 
temperatures by more than 10°C, and, thus, the air can be 
unstable. This atmospheric instability increases the air-sea 
thermal exchange and results in the high salinity layer of 
the sea becoming unstable. Hence, convective sinking occurs 
and the thermocline disappears. 

This process, which 1s particularly noticeable in 
the Rhodes-Cyprus area, results in the formation of a 
Levantine Intermediate Water (LIW). It is a mixture of over- 
turned surface water, Atlantic Water (which was sheltered 
by the summer thermocline) and the water formed in the 
previous year. This homogeneous water body is usually found 
initially at a depth of 150m, near its source, down to 400m. 
Following a cyclonic gyre in the Levantine Basin, the LIW 
mews through the Strait of Sicily, continuing towards the 
Western BaSin and into the Atlantic Ocean. Between its 
Bemece and the Strait of Sicily, it flows isentropically, 


with a change of its characteristic temperature and salinity 
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mmeimelcre © 39. 5a at ats Source to 15.8°C, 38.90% in 
m@emronwan Sea and 14.0°C, 38.74% in the Strait of Sicily 
fezencurt, 1976). 

As described above in the Bathymetry Section, the 
Maltese Sill provides the most obvious passage for this 
meemen tlow, which has been estimated to be 10 million m’s7', 
with an excess of 4.5% of in-flow over out-flow. Miller 
[1972] suggests that the westward flow of LIW is more prom- 
inent in the northern part of the Strait, i.e., between Malta 
and Sicily, while the eastward flow of Atlantic Water 
dominates the southern part. 

In addition to this main flow pattern, there is some 
evidence of summertime exchange between the Ionian Sea and 
the Adriatic Sea [Zore-Armanda, 1969]. Due to the continental 
Character of the narrow, shallow Adriatic Sea, its water is 
lighter than Ionian Sea Water, with a noticeable input of 
fresh water from the northern rivers of Italy. This surfi- 
cial water flows south into the Ionian Sea, and a compensation 
current is formed in the intermediate layer, with Mediterra- 
nean Water flowing into the Adriatic Sea at around 300m 
depth. 

meee Frontal structures in the Ionian Sea 

The Western Ionian Sea is an area Of Mixing. Within 
the upper 100m of the water column, which includes the 
thermocline, Atlantic Water spreads over the area while LIW 


Flowing underneath is forced by topographical features 
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maeougia the narrow Stll of Malta (Fig. 3). As can be 


expected, thermal and salinity fronts occur as quasi- 
permanent features, although marked by great variability in 
position and strength. Several studies of these fronts have 
been made [Woods, 1972; Johannessen et al., 1971], some of 
the features noted are described below. 
In May, July, and August, a front (called the Maltese 
Front) was found along the continental slope east of Malta, 
extending about 130 km from the southern tip of Sicily. 
Cross-front horizontal gradients reached 2°C over 33 km 
(20 nmi), with extensive vertical shear over a depth of 100m 
due to up- and down-welling processes along the interface. 
Surface water flowed through the front into the Ionian Sea. 
The most characteristic features of the Front seem 
to be: 
its position is associated with bottom topography 
a Salinity and temperature minimum lies to the west of 
the salinity and temperature front, and 
temperature inversions are highly probable [Briscoe 
Se ale, 1972). 
Two studies of the Ionian Sea made in the late 
Summer and fall [Levine and White, 1972; R. R. Miller, 1972] 
Showed the existence of an oceanic, 1.e., deep water, frontal 
structure, located several tens of kilometers east of the 
Maltese Front. Its main discontinuity was manifested in the 


presence of warm anomalies of LIW extending to a depth of 
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400m and presumably moving westward. The frontal features 
did not always appear in the sea surface fields, but they 
were mainly noticeable at depths between 50 and 400m. 

It was later suggested that the Maltese Front under- 
goes seasonal variations and propagates eastward, from summer 
to winter, into the Ionian Sea [Johannessen, 1972]. 

Woods et al. [1977] have discussed a model of verti-~ 
cal circulation for fronts in this area which displayed the 
upward and downward displacements of water tongues. The 
shape of isotherms was described and related to features 
found commonly in the temperature transects, depending on 
the orientation of the cross-section with respect to the 
meental axis. 

Woods [1977] also examined the analysis made by 
Briscoe et al. [1974] to show that the assumed linear frontal 
zone, interpreted by Briscoe et al. as extending over 100 km, 
might possibly have derived from an erroneous interpretation. 
Woods emphasized that a 50cm/s southward drift was not taken 
into account, and that what was described as sections of an 
elongated front could just as well have been the repeated 
crossing of a water mass boundary of limited extent. As an 
alternative, Woods showed a preference for the concept of 
widely spread, but small scale, frontal structures, imbedded 


in a turbulent upper ocean. 
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A. BATHYTHERMOGRAPH DATA 

The present study was initiated when a set of one hundred 
and sixty seven expendable bathythermographs (XBTs) were made 
available by FNOC after completion of the US Navy exercise, 
SHAREM 38, in the Western Ionian Sea in the middle of June, 
mo 30. 

The first step in evaluating this data set was to plot 
each temperature profile, its geographical location, and to 
contour temperatures at various depths. 

1. Vertical temperature profiles 

Vertical temperature profiles provided a visual 
estimation of the quality of the 167 XBTs. The quality of 
each XBT was examined according to the following 
Seiaracteristics: 

the number of points in each profile and maximum depth 
of the sounding. 

- the smoothness of the profile, 1.e., did the sounding 
include the small scale features generally expected, or 
had it been excessively "smoothed"? 
the presence of discontinuities, 1.e., did they repre- 
sent an actual thermal feature or merely a defective 


sounding? 
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The data set was then separated into three categories. The 
best quality category included those profiles which extended 
to at least 150m, had a large number of depth/temperature 
pairs and resolved reasonably well the 30 to 100m layer. 
Rejection of profiles was made on the basis of limited verti- 
cal extent (less than 40m), or an obvious anomaly. This 
constituted the category of rejected profiles. The remaining 
profiles constituted the second quality category. The best 
quality, second quality and rejected profile categories 
contained 74, 70 and 23 profiles, respectively. Figure 4 
shows two XBT profiles that were classified, in spite of 
their similar shape, as best and second quality due to the 
absence, in one of them, of the typical inversion feature. 

No effort was made to differentiate surface-launched 
SXBTs from alr-launched AXBTs and the appelation XBT was 
used for both types of bathythermogram. 

2. Geographical distribution 

The geographical distribution of XBTs was found to 
be concentrated in two well-defined areas. This spatial 
pattern was supposedly due to a change of operational area 
during the exercise. Zone I was defined as covering 33°30' 
meme N, 17° to 19° E. Of the 67 soundings from Zone I, 
45 were first quality, and 19 were second quality soundings. 
Meme ii was defined as 35° to 38° N, 16°30' to 19° E, thus 
GOvering four times as much area as Zone I. Of the 85 


soundings from Zone II, 23 were first quality and 47 were 
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second quality soundings. Finally, 15 soundings were widely 
spread over the rest of Western Ionian Sea. Table I summari- 


zes these facts. 


TABLE I 


Cleltievyeana geographical distribution of XBTs 


Best Second rcp lott— 
quality equality Rejected fecal able 
zone I 45 ng 5 67 64 
Zone II 23 47 iS 85 70 
Other 6 4 a 15 A) 
rotal 74 70 23 167 144 


Some disproportion seems to arise regarding the 
memcent of first guality profiles in Zones I and II, 
respectively 67% and 27% of the total number. This dispro- 
PeeelOon is certainly due in part to the inherent subjectivity 
@meene Quality classification, but it points out that pro- 
files in Zone II, which were more frequently classified as 
“exceedingly smoothed out, hence second quality" might in 
reality be smoother. Hence, Zone I could be expected to be 
an area with more variability in the vertical structure: 
inversions, sheets, layers, etc. In spite of this subject- 
ivity, the quality classification was retained as a valuable 


mein thermal field contouring. 
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pee lame ALStrtourron 
The 70 XBTs in Zone II were acquired between June 6 
and 24, with some concentration (64) from 13 to 20 June. 
Similarly, the 64 KXBTs in Zone I were acquired between 
June 6 and 15, with some concentration (55) from 8 to 12 
wune. 
The two data subsets: 
Zone I - 8 to 12 June 
zone II - 13 to 20 June 
were finally selected for independent. thermal analysis based 
eae ehe following criteria: 
- good apparent quality of the individual profiles 
- relatively high spatial density, and 
- potential use as a synoptic set as discussed below. 
Figure 5 summarizes the space and time distribution of 
profiles for the two selected subsets. The use of ‘Zone I' 
and 'Zone II' in what follows 1s indicative of both the 
memeoral and the spatial limits defined above. 
4. Thermal analysis 
An analysis was made for temperature structure at 
the surface, 100, 200, 300, and 400m. The thermocline depth 
and inversion depth were mapped. Confidence was given to 
Singular points according to the gualitative estimation of 
each XBT, i.e., a best or second quality profile. For 


example, the absence of an inversion in a second-quality 
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sample was not considered significant when neighbouring 
first-quality soundings showed such a feature. 

A cross section was evaluated to a depth of 400m 
where thermal features presented the greatest horizontal 
variation. A profile of thermocline depth and inversion 
depth was superimposed on this cross-section. 

Discussion of this analysis, its validity and the 


results are presented in Chapter VI and VII. 


eee CORRELATION COMPUTATIONS 

A statistical approach was used in both Zones I and II 
to determine the horizontal structure of the perturbation 
field at different levels. It included the following 
steps: 

1. A perturbation field was generated as the difference 
between a given XBT profile and the mean profile for 
the area uSing a vertical resolution of 5m. 

2. XBTs were then taken 2 by 2, and a cross correlation 
coefficient was computed over three different layers: 
0 to 60m, 60 to 150m, 150 to 400m. These layers were 
selected as representative of different physical 
processes in three regimes: the well-mixed layer, 
the zone of inversions, and the quasi-isothermal layer, 
respectively. 

3. The cross-correlation functions were evaluated as a 


function of the spatial and temporal interval between 
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the two profiles, in an attempt to validate the 
assumption of synopticity. 
4. The cross-correlation functions were also analysed as 
a horizontal field with respect to latitudinal and 
longitudinal spacing between the XBT pairs. 
Details of the computational technique are given in 


Appendix A. 


ee SATELLITE IMAGES 

A set of daily Advanced Very High Resolution Radiometer 
(AVHRR) satellite images was obtained from the French 
Meteorological Office, covering the area of interest from 
mid-May to the end of June. The quality and cloud coverage 
varied widely throughout this period and ultimately 14 
pictures were selected which gave a good estimation of skin 
SST gradients. These pictures were sometimes used in pairs, 
when they were not too separated in time and showed no 
evolution of the characteristic features, to make a comvosite 
image. The resolution between two tones of gray was 0.5°C 
and no correction was attempted for atmospheric absorption. 
Hence, the contours were taken as representative of SST 
contours, whenever atmospheric effects did not apparently 
fie the picture. 

Since the satellite images obtained were direct enhance- 
ments, for oceanographic purposes, of raw IR pictures, they 


Zee nad a grid distortion corresponding to the particular 
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geometry of the satellite (NOAA-6 or TIROS-N) scanning 
procedure. A computer program was written to rectify each 
immage on a geographical grid. The numerical technique 
mie luded : 
ie Gigitization of the image contour, 
2. initialization of the scanning grid for each 
Barticular orbit, 
3. point by point inversion using an 8-noded isopara- 
metric mapping, and 
Meee Output and plotting. 
All of these steps were executed on the Tektronix 4081 of 
the Mechanical Engineering Department. The restritution was 
believed to have a sufficient accuracy to locate thermal 


features within a few kilometers. 


D. WIND FIELD ANALYSIS 

A set of daily wind observations at 1200Z (from 1-20 
June) was obtained from FNOC, with wind vectors evenly 
@motributed at 16 grid points throughout the study area. 
Surface pressure analyses with the same periodicity were 
received from the French Meteorological Office. 

The wind stress at the surface was computed for each 
grid point and each day as 7 = p.Cpulul, = being the 
density of the air, u the wind Wellectey sUsindel+5 20" fOr 
the dimensionless drag coefficient C,. The curl of the wind 


stress was computed for 9 intermediate grid points, using 


34 





an elementary finite difference scheme. Using a result of 
Ekman theory based on the conservation of water mass and an 
evaluation of the mean horizontal mass transport, the verti- 
cal velocity at the bottom of the Ekman layer was computed 


as: 


1 = 
— (¥ x Tt) 
ote 


= 
il 


where 0 is the density of sea water, f the Coriolis parameter 


and V x the rotational operator. 
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IV. WEATHER SITUATION 


me DESCRIPTION 

During early June, 1980 the synoptic weather situation 
was mostly determined by the Azores anticyclone, extending 
over western Europe, and a series of depressions over the 
Balkans and central Europe. Winds over the Ionian Sea were 
eastward, 5 to 10m/s (Fig. 6a). This situation remained 
steady until June 6, when the Azores high regressed anda 
series of depressions developed over Spain, the British 
Isles and Algeria. By June 8, this latter center of low 
pressure had caused the wind to veer to northward, 5 to 7m/s, 
and with a speed increase to 10m/s on June 10. During that 
two-day period, some cloudiness was present over the 
Ionian Sea. 

After June 10, an anticyclone started developing over the 
central Mediterranean Sea, and winds were moderate and un- 
stable in direction through June 12. On June 13 and 14, 
winds were steadily westward, 10m/s, while a strong depres- 
Sion developed over the Bay of Biscay and passed over 
northern Europe. The following period, June 16-20, was more 
typical of the seasonal pattern, with a moderate wind regime 
over the Western Ionian Sea. Winds were south- to south- 
eastward, 2m/s on June 17, increasing to 8m/s by June 20. 
Air temperatures, aS reported by the Malta weather station, 


were over 24°C throughout most of the month. 
36 





eee LMPACT OF WIND STRESS 

Dynamical theory indicates that the surface wind stress 
increases turbulence in the upper layer and contributes to 
a downward redistribution of heat, i.e., deepening and cool- 
ing of the mixed layer, or steepening of the thermocline. 
Moreover, transient winds can generate inertial oscillations 
[Pollard, 1970], an effect which has been observed in the 
Mediterranean Sea [Perkins, 1972; Mittelstaedt, 1973]. 

This is most noticeable during the period of calm following 
a period of high winds [Gonella, 1971]. 

ie hkman “pumping” or “Suction” due to the curl of the 
wind stress was evaluated (Fig. 6b) as a time series of 
daily vertical velocities, computed at the bottom of the 
Ekman layer (see Chapter III). The most striking event in 
this series is a very strong upwelling induced during the 
period 10 to 12 June, with mean values exceeding 0.02 cm/s 
(l5m/day), a very high value indeed. This was followed by 
a more moderate, but still noticeable, downwelling during 
mee period 13 to 15 June. 

It is not known, from available observations, whether 
maes vertical "pumping" or "Suction" acutally interacted 
with mixed-layer processes or not. From theoretical consid- 
erations, it would raise or depress the water column, hence 
the thermocline, without heat redistribution. A very strong 
Moewelling, possibly coupled with strong wind-mixing, might 


rupture the thermocline and account for the presence of 
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intermediate water at the surface. Since these estimations 
were not applied as corrections to the XBTs, any vertical 
displacement of isotherms by such an Ekman "pumping" or 


Memection" contributed to increase the noise level in the 


data. 
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V. SPACE AND TIME VARIABILITY 


To validate the analysis of the thermal field, it was 
necessary for all XBTs taken in each of the two periods to 
be established as essentially synoptic with respect toa 
temporal scale of 8 days and a spatial scale of about 30 km. 
Eight days was the duration of the sampling effort in Zone 
II, the longest of the two data subsets, and 30 km was an 
average sampling interval in the densest area. Among the 
mechanisms which could possibly distort the analysis, l.e., 
falsify the synopticity assumption, were horizontal and 
vertical advection due to mean flows, internal waves, etc., 
surface heating, interleaving and mixing of water masses, 
all of which might have occurred during the sampling interval. 
It will be shown in the following sections that the above 
mechanisms did not violate the assumption of synopticity at 


the temporal and spatial scales specified. 


A. VARIABILITY OF SURFACE FEATURES 
1. Characteristics of satellite images 

The series of seven satellite composite pictures 
(Fig. 7-13) was found to have the following features in 
common. 

A persistent tongue of cold water oriented eastward 
from the sill, at 35°20'N, was distinguishable and observed 
as far east as 18°E. It presented a varying shape, either 
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Pensa gnem~¢mig. Ll and 13) or bent southward around longitude 
17°20'E (Pig. 7, 10, and 12). In some cases it turned into 
a clockwise circular pattern, suggestive of an anticyclonic 
pee (Figs. 8 and 9, and possibly Fig. 12). 

The other dominant feature was a tongue of cooler 
water issuing also from the Maltese Sill, which tended to 
bend southward along the Medina Bank. It then followed the 
Shelf break, along 16°E and bent eastward again along 34°N. 
This feature was seen on all seven satellite pictures and 
was separated from the tongue of cooler water mentioned above 
by an area of warm water, sometimes in anticyclonic eddy- 
type shape, centered approximately at 35°N, 17°E (Figs. 8, 
ever). 

In the northern part of the Ionian Sea, a persistent 
tongue of cooler water was observed flowing southward from 
the Strait of Messina. Upwelling structures could be seen 
along the eastern coast of Sicily, with occasional extensions 
offshore (Fig. 11). Within Zone II, the surface features 
seemed to be quite stable, and surface temperatures were 
notably warmer than in the southern part. 

Based on the composite pattern from 4 to 30 June 
(Fig. 14), there are commonly two branches in the sea-surface 
temperature field, as found in earlier Airborne Radiometer 
Temperature (ART) studies, e.g., Briscoe et al. [1974; 

Fig. l]. Using the same ART technique, Saunders [1972] 


established the existence of a tongue of cold water which 
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corresponded exactly, in season and geographical location, 


to the northern most of the two cold tongues found in the 
present study. Finally the branching of the cold tongues 
can be traced also from a salinity analysis [Morel, 1972; 
meee S| along 15°E, 1.e€., across the sill, in which two 
shallow salinity minima (Atlantic Water) are found 45 km 
apart. 

All these factors suggest a relation between this 
typical shape of SST contours and the circulation of water 
mmaeough the Strait of Sicily. 

2. Surface temperature and water mass circulation 

As noted by Saunders [1972,p473], "there are clear 
connections between the surface circulation and the surface 
temperature field." Although it was not possible to cali- 
brate the satellite imagery, it was clear that such patterns 
of cooler and warm water, with their large scale and persist- 
ence, were mainly due to advection and not to loss of heat 
@emreaistribution of heat in the ocean. This observation 
was used by Miller [1972, pl4] to assert that 

transient as the surface temperature may be, advective 
Se@rculation could orient the distribution of surface 
temperatures in such a way as to indicate the super- 
position of intruding water masses from the anomalous 
features of the distribution, the apparent organization 
and horizontal coherence of isotherms, and the rapid 
changes of temperature with short distances. 

It is also noted that double-tongued flows are not 


atypical downstream of the Straits of Sicily and Gibraltar, 


and, in a more general frame, of vortex patterns. 
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The SST pattern was thus interpreted as two flows of 
cold water from the Strait of Sicily, the southern tongue 
following a topographically controlled path along the shelf 
edge, while the northern tongue curled into an anticyclonic 
im LOwW. 

The eddy structure which existed on June 14 and 15 
(Figs. 8 and 9) did not appear on June 22 and 23 (Figs. 10 
geet!) , but was found again on June 24-25 (Fig. 12). Hence, 
this pattern can be apparent, disappear, and appear again 
in 4 to 8 days. This time scale is typical of a nominal 
weather cycle; however, no definitive study has been made 
of this time scale of variability in the Ionian Sea. As an 
explanation, it could be hypothesized that the summertime 
heating of the upper ocean masks its deeper structure, until 
a weather event causes sufficient mixing to unveil this 
thermal structure. Hence, the atmosphere would not neces- 
Sarily provide the generating force of these oceanic thermal 
features, but would act as a revelator. 

eee fhe etfect of advection 

The velocity of water masses in the area has been 
estimated as follows [Saunders, 1972] from the MILOC-68 data 
Sec : 

Emonguco’N, 1.€., in the tongue of Atlantic Water, 
currents are eastward at 10 to 20 cm/s 
Pmouna 37°N, 18°E, Ehe Currents are weak (5 to 8 cm/s) 


and northeastward. 
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These current estimates were based on geostrophic velocities 
computed relative to the 200dbar level. No indication of 
westward flow was found at the level of LIW. The surface 
flow seems to be rather sluggish in most of Zone II, while 
the mean flow is stronger south of 35°30'N. In the vicinity 
of frontal features, velocities are typically greater, with 

a magnitude of the order of 30 to 50 cm/s. Such velocities 
would occur in high shear situations, along water mass 
boundaries, which in this region are spread north-south, as 
shown by Johannessen [1971; 1977] and others. Givena 
spatial scale of 30 km and typical advection speeds of 5 cm/s 
(general case) and 35 cm/s (frontal areas), a tentative 
Synoptic time scale was computed as the ratio of length scale 
to surface current speeds, yielding values of, respectively, 
seven and one day. 

From this evaluation of the effects of advection only, 
it was inferred that the XBT data could be considered as 
Synoptic over a seven-day interval as long as they did not 
Show the characteristics of a fully formed frontal structure. 
In the presence of such a feature, a one-day interval should 
be held as a limit of synopticity, unless some corrections 
are made to account for advection, due to the high velocities 


Ememene variability of mesoscale features. 


B. EVALUATION OF TEMPORAL VARIABILITY 
This section describes an attempt to evaluate the import- 


ance of temporal variability in each of the two data subsets. 
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A mean temperature profile was computed for Zone I, as if all 
XBT profiles were co-located. This average profile is refer- 
red to, later in the text, as the areal mean profile. 
Similarly, a mean temperature profile was computed for each 
day, and the five daily mean profiles from the June 8-12 
period were evaluated as a time series. Since it was diffi- 
cult to separate visually such profiles when plotted on top 
of each other, the perturbation profiles, i.e., the arith- 
metic difference between each daily profile and the areal 
mean profile, were derived. The analogous procedure was 
Bemeormed in Zone II (Fig. 16). 

In Zone I, there was a marked difference in the 8 June 
mean perturbation profile: it was significantly cooler than 
the average in the well-mixed layer, down to 100m. This was 
believed to express a geographical characteristic, 1.e., a 
local cold area which could be the southernmost part of the 
tongue of cold water seen in Figure 7. It 1s not considered 
as representative of a temporal phase in the area. This 
interpretation is validated by the observation that the mean 
Meederon (Fig. 5) for June 8 is Significantly different from 
those of the other four days. 

The mean profiles for June 9 and 10 are quite similar, 
with no major deviation from the areal mean profile. 
memrmrarly, the profiles for June 11 and 12 are within 0.25°C 
of the areal mean profile and can be regarded as a pair of 


Similar data sets. But a comparison of the two vdairs, 
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June 9 and 10 with June 1l and 12, in the upper 25m, shows 
that a Significant (0.5°C) average heating occurred after 
June 10. This can be well related to the weather pattern 
described in Chapter IV, where June 10 was noted as a 
separation date between a period of strong winds with cloudy 
conditions, and a calm period with intense heating. In 
addition, the well mixed layer between 25 and 75m revealed 
markedly variable temperatures, as compared with the areal 
mean, although with no obvious time dependence. From this, 
it could be expected that spatial analysis of the thermal 
field in this upper layer would be difficult. 

Temporal variability in Zone II was difficult to assess, 
because the initial assumption that all XBTs were co-located 
was far from true. On the average, the displacement of the 
sampling ships from June 13 to 20 was northward, about 30 
km/day, although erratic. This 1s referred to below as the 
mean sampling speed. There was an apparent heating of the 
surface layer, with mean surface temperatures increasing in 
eeemtOollowing time sequence: 13, 14, 16, 15, 17, 20, 19, 

18 June (Fig. 16). At this point, it was not possible to 
separate spatial from temporal effects without an alternative 
source of data. The satellite images could not help, since 
they were not calibrated. Moreover, as the spatial density 
of the XBTs was low, no small scale resolution could be 
expected in Zone II, in contrast with Zone I. Hence, it was 


not possible to separate the real heat gain in the surficial 


45 





layers, which occurred most probably, from the apparent 
warming due to the mean direction of sampling, from cooler 
to warmer waters. In the subsurface thermal field, some 
grouping of the perturbation profiles is observed around 
200m, with June 17, 18, and 20 being warmer than the areal 
mean and June 13, 14, 15, 16 and 19 being cooler. Given the 
depth of this grouping, this is not believed to represent 

an effect of surficial heating, but more probably an indica- 
tion of the spatial distribution of water masses or some 


internal process, i.e., large-scale wave or eddy. 


oe)6 LEMPORAL VERSUS SPATIAL CORRELATION 

The following section describes an attempt to use the 
cross~correlation computations (see Chapter III and Appendix 
A) to enlighten the respective roles of temporal and spatial 
variability. The following reference system is defined: 
one axis describes the time difference, in increments of 
days, between XBT pairs; the other axis describes the dis- 
tance in increments of 10 km between pairs. In a perfect 
steady-state situation, the cross-correlation coefficients 
between XBTs would vary only with distance. Hence a plot of 
these coefficients in the reference system would show 1so- 
Beerelation contours lying parallel to the time axis (Fig. 
17a). Conversely, the variation with time of a homogeneous 
temperature field would be revealed by iso-correlation 


Bemeours lying parallel to the distance axis (Fig. 17b). 
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1. Observation in Zones I and II 

In Zone I the correlation decreased with distance 
(Fig. l17c) but it was quite constant with time. This is 
dramatically illustrated by the fact that two profiles 
separated by 30 km correlated at a 0.5 level, even when 
sampled at a four-day interval. 

By contrast, in Zone II there was a clear dependence 
of iso-correlation contours in both time and space (Fig. 17d). 
The 0.3 correlation level existed as far as 80 km for a pair 
of XBTs acquired on the same day, and as far as 30 km for 
XBTS acquired two days apart, but not for any larger time 
Semaration. Similarly, the zero-correlation contour, which 
occurred at a distance of 150 km for same-day samples, was 
found at a four-day separation for co-located samples. 

From Figure l7c Zone I was essentially in steady 
state over the five-day period, even in the presence of 
strong advection. This concept of synopticity applies 
relative to a length scale of, typically, 30 km. It was 
consistent with satellite pictures which showed an average 
temporal scale of three-to-eight days for the tcongue of 
Atlantic Water and related patterns. 

In Zone II, as was seen in the preceding section, 
it was not possible to separate the various effects which 
may have affected the temporal vs. spatial correlation: 
advection, heating processes and mean speed of sampling. 


From Figure 17d, the speed scale is 150 km/4 days for the 





mee correlation contour, 1.e., 37 km/day, and 90 km/3 days 
for the 0.3 correlation contour, or 30 km/day. These two 
values are fairly close to the mean sampling speed. It is 
believed that the temporal distribution of samples over an 
eight-day period in Zone II represents the main cause in the 
loss of correlation with time. Hence, the physical processes 
involved could be expected to have a longer time scale, 
implying a true synopticity in Zone II, even though it should 
ideally be analysed as a time succession of smaller synoptic 


"analysis spaces." 
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VI. THERMAL ANALYSIS IN ZONE I 


meee DEoCRIPTION OF THERMAL STRUCTURE 

This section describes the thermal structure found in 
Zone I and refers to Figs. 18 through 24. 

There was a core (diameter of ca. 30km) of warm surface 
water (20.2°C) centered at 34°20'N, 18°15'E surrounded by 
water approximately 0.5°C colder (Fig. 18). The northwestern 
part of the area was somewhat colder (19.2°C), and seemed to 
correspond to the cold water tongue, or its distant influences, 
intruding eastward. The distribution of XBTs in that part of 
the zone was not adequate to resolve the two tongues of water 
appearing in the satellite imagery. 

The SST analysis was made much easier by using an "“inter- 
Memacea Sot’: the Expanded Ocean Thermal Structure (EOTS) 
model run by FNOC provides a grid analysis of thermal struc- 
ture from input XBTs, following a scheme that is briefly 
described below (Chap. IX). A by-product of this analysis is 
a geographical interpolation of the SST analysis, which 
provides an analysed SST at each XBT location interpolated 
from SST values analyzed at EOTS grid points. This inter- 
polated SST is generally believed to smooth out small scale 
Structure and only retain well correlated features. In this 
respect, it seems to overestimate the warm water core 


described above rather than average it out (Fig. 19). 
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The analysis at 25 and 50m did not show any consistent 
pattern. Since these depths are located above or in the 
average thermocline, it is not surprising to find extreme 
Variability there. 

Below the thermocline, the analyses at 200, 300, and 
400m showed a warm core area 35 km eastward of the equivalent 
surface feature; this displacement increased with depth 
(Pigs. 18-22). Average horizontal gradients through this 
SescOntinuity were of the order of 1°C/50 km. 

The thermocline depth, defined as the depth with a maxi- 
mem temperature gradient, has a zone of transition between 
a shallow depth (less than 40m) in the western partof the 
area to a deeper depth (greater than 60m) in the eastern 
part. The thermocline depth had strong variations in the 
region of the warm anomaly, including an area where the mixed 
layer seems to vanish. 

Most XBT profiles in the best quality subset of data 
showed inversions, i.e., the presence of warm water under- 
lying cool water. Such a feature is typical of regions where 
the salinity gradient has a dominant role in the density 
Stratification. The inversion depth increased from the 
northwest corner (50m) to the southeast corner (125m) (Fig. 
23). Again, the warm anomaly was associated with perturba- 
tions in the inversion depth, which decreased sharply over a 
Short horizontal distance (from 100 to 50m in about 10 km), 
just east of the warm anomaly and precisely where the thermo- 


meme seemed to vanish. 
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PEt aniseemewas enstruced from ten XBTs spanning 120 kn, 
with a greater concentration of XBTs in the perturbed area. 
The characteristics of the horizontal and vertical tempera- 
ture fields are summarized below (Fig. 24): 

- a deep (200m and more) warm anomaly was present, which 
was characterized by a large isothermal (15°C) body of 
water, whereas the 15°C isotherm is found above 150m 
depth throughout the remainder of the area (Fig. 25). 

*" a strong thermocline was present west of the warm 
anomaly; 

* the temperature variability was relatively great in the 
surface layer; 

" jlayered inversion zones were present; 

* relatively warm surface water was present on the western 


edge of the warm anomaly. 


Pee COMPARISON OF SATELLITE AND XBT SEA-SURFACE ANALYSES 

An attempt was made to compare thermal features observed 
in the satellite IR images (Figs. 7 and 8) with those seen 
in the surface XBT analysis. The comparison did not prove 
very successful, due to a paucity of XBTs where the satellite 
imagery presented the most characteristic eddv-shaped feature. 
However, the satellite IR image (Fig. 8) shows a cool (1.5°C) 
anomaly? where the interpolated XBT SST shows a warm anomaly. 


taithough this satellite image 1S a composite made of two 


actual images, onlv the June 14 picture showed that 
particular feature. 
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Several factors suggest that the satellite anomaly might be 
due to some atmospheric effect: its size, the fact that it 
appears on one IR image only (no persistence), and its 
transparence to some of the surrounding oceanic features. 
The development of marine fog or stratus from a warm water 
patch has been observed and modeled [Rogers and Hanley, 
1980], and could explain the phenomenon described above. No 
attempt was made to verify this hypothesis from local 


atmospheric observations. 


fo CORRELATION COMPUTATION 

It can be recalled from Chapter III that a correlation 
coefficient was defined to allow for quantitative comparisons 
of a pair of XBT profiles. The number-density contours of 
XBT correlation pairs in each 10 by 10 km square generally 
have a E-W stretching (Fig. 26). Hence more profile pairs 
had a east-west than a north-south orientation. 

The cross-correlation of near-surface temperatures 
Memo. 27) haS a wavy pattern of significant (+0.5) positive 
correlation in the N-S direction, at distances of 25 and 60 
Km. Zero-crossings occurred at 30 km in the E-W direction 
and 40 km in the N-S direction. 

The N-S pattern exhibited in the surface layer was even 
more notable in the intermediate (60-150m) level (Fig. 28). 
but there was also an area of significant negative correla- 


tion at a distance of 70 km east and west or the center. 





Similar patterns were found in the deep (150-400m) layer 
more, 29). 

Due both to measurement noise and the impossibility of 
resolving details within the averaging grid size of 10 km, 
all three correlation diagrams had values less than +1.0 at 
Zero Spatial lag. 

The thermal structure in the area was clearly anisotropic. 
Choosing the deep layer as the most stationary, two distinct 
patterns are distinguished: 

- in the E-W direction, where the zero-crossing yeilded 
a correlation scale of 50 km, physical phenomena seemed 
to occur with a typical separation distance of 200 km. 
In fact, the correlation fields (Fig. 27) displayed in 
a different fashion the most apparent features described 
in the thermal analysis section, in which one could 
observe (Fig. 18): 
- a warmer area in the east 
- a colder area in the west 
- a strong thermal gradient area with lsotherms 

Grilented in the N-S direction. 

meen the N-S direction, the spatial scale of variability 
was larger. With a Separation as large as 80 km, XBTs 
still showed a Significant positive correlation. 
Superimposed on this positive correlation field was a 
wave pattern, with a length scale of approximately 40 km, 


1.€., Similar to the length scale in the E-W direction. 





Hence, similar wave or eddy processes occurred in all 
directions, though they were superimposed on a strong 
positive E-W pattern. This result is also in good 
general accordance with the thermal structure analysis. 
A natural axis of symmetry, along 340°T, 1S apparent 
and could represent, with the general E-W axis, a 
preferential axis for advection or baroclinic wave 


propagation. 
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VII. THERMAL ANALYSIS IN ZONE II 


The location and the circulation of water masses in Zone 
II are basically similar to those in Zone I, as described in 
Chapter II. Zone II's bottom topography is somewhat flatter 
than Zone I's. Zone II is more sheltered from the Western 
Mediterranean Sea, but it seems to be more exposed to the 
westward LIW flow, as well as possible exchange with the 
Adriatic Sea. The data density was one profile per 1650 sq 
km, aS compared to one per 550 sq km in Zone I. These 
factors led to a separate analysis in Zone II, while expect- 


ing a lesser resolution than in Zone I. 


Swe ORoCRIPTION OF THERMAL STROCTURE 

The SST pattern demonstrates a general NW-SE orientation 
of isotherms, with temperature increasing in the NE direction 
Mero. 30). There was a “tongue” of cool (1°C) water along 
Peemewnich bent southward along 17°30'E. It could not be 
accurately matched with satellite IR images, although its 
sourthern boundary correlated quite well with the eddy-shaped 
Structure of Figs. 8 and 9. The XBT data were sparse SE of 
this feature; hence, the positive temperature gradient is 
mimey Significant to the NE. 

The interpolated (EOTS) SST field (Fig. 31) shows the 


Same basic features, and, in particular, the 20° and 21°C 
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isotherms are almost identical to the XBT-derived SST. 
However, the EOTS SST does not accurately depict the 22°C 
isotherm, nor does it depict well the cool anomaly near 38°N, 
17°40'E in the north of Zone II. Analyses were conducted at 
the 100, 200, 300, and 400m levels and showed little change 
with depth. Hence, the subsurface temperature field is only 
shown here at the 100m level (Fig. 32). 

At all levels, there was a warm anomaly in the NE of Zone 
II, with a general orientation of isotherms more N-S than in 
the surface analysis. There was a band of cool water at all 
levels between 16°E and 17°E, and it seemed to widen as it 
reached 36°N. Here, it joined the cool area centered at 
35°30'N, 17°30'E, which was the subsurface expression of the 
surficial tongue of cool water described above. This cool 
tongue of water could be interpreted as a local frontal 
structure, whereas the elongated band extending northward 
could indicate exchange with the Adriatic Sea, at an inter- 
mediate level. 

The correlation technique described previously was applied 
to Zone II. The horizontal cross-correlation diagrams were 
of little value, due to the high "noise" level in the data. 
The analysis was, hence, made non-directionally (Fig. 33), as 
an averaging of the time vs. space correlation (Fig. 17d) along 
the zero lag time axis for the eight-day period. The zero- 
Seessing correlation distance was 50 km in the 150 to 400m 


layer, in good agreement with the results in Zone I. 
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Veins toon ACOUSTIC FEATURES 


The goal of the acoustical analysis was twofold. First, 
some particular oceanic features might be best revealed by 
sound speed profiles. This is usually true when looking at 
extrema; e.g. sound speed maxima or minima, maximum gradi- 
ents, etc. Second, some understanding of the impact of these 
features on ASW detection capabilities might develop; e.g., 
how tactical parameters (surface duct, deep sound channel, 
critical depth, etc.) relate to easily accessible oceanic 


varlables. 


A. PRELIMINARY REMARKS, DEFINITIONS 

The spatial and temporal distribution of XBTs has been 
described extensively above. The following analysis does 
not deal explicitly with physical processes, and the data 
are presented regardless of previously defined time intervals, 
Subzones, etc. Moreover, some previously ignored XBTs, 
Sparsely distributed southwest of 36°N, 17°E, have been 
incorporated. 

The Leroy [1969] equation was used to convert the XBT 
orofiles into sound speed profiles (SSP). Because of the 
absence of salinity data associated with the temperature 
profiles, a single climatological salinity profile from FNOC 


archives was used. The effects of this approximation were 
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checked and it proved to have no significant effect on the 
shape of the SSP, nor on the depths of its characteristic 
features. The largest deviation in absolute sound speed 
from a test profile was 1 m/s at the surface and deviations 
decreased with depth. 

The more typical sound speed profile (SSP) characteristics 
were defined below (Fig. 34) in terms of their nature (gradi- 
ent, extrema), their physical origin (inversion, mixed layer), 
or their effect on sound propagation (duct, sound channel): 

- a near-surface positive gradient” which 1S due to a 
surficial isothermal layer and relates to the acoustical 
surface duct. 

* a strong negative gradient, which occurs where the 
temperature decreases most strongly with depth, and 
defines the thermocline depth. 

- the main sound channel, or summertime channel, whose 
existence 1s due to the seasonal heating of the surface 
layer above a deep near-isothermal water mass. Its 
axis cOincides with the absolute sound speed minimum. 
Its thickness is the difference between the depth of 
the absolute near-surface sound speed maximum and the 
depth at which this same sound speed is encountered 
below the axis. 


2 : 
the sign convention adopted here 1S common to underwater 


acousticians, (as opposed to the oceanographers' convention), 
whereby depth is positive downward. 
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one or several secondary sound channels, which are due 

to the presence of thermal inversions or isothermal 

water masses, and are characterized by an axis, a thick- 
ness and a mean gradient. Axis and thickness are defined 
as above, but the sound speed maximum is, in this case, 
the nearest (to the minimum) relative sound speed 
maximum. The mean gradient is defined as 2AS/AZ, where 
AS is the difference between local sound speed extrema, 
and AZ is the channel thickness. 

a deep positive gradient due to the dominant effect of 


pressure on sound speed. 


pee VLSTRIBUTION OF VARIABLES 

Subareas were re-defined to increase the contrast from 
one subarea to another, in order to facilitate the study of 
the distribution of sound channels. The water mass properties 
and the topography were taken into account. This led to the 
definition of five subareas, A through E, A being located on 
the continental shelf, B and C in Zone I, D in the coastal 
Part of Zone II and E in the oceanic part of Zone II. 

The thermocline depth was a systematic variable, l.e., 
its depth could be defined almost unambiguously and changed 
meeciiy from one location to another (Fig. 35). Its mean 
value was about 40m in Zone II, and varied between 30 and 
70m in Zone I. Greater variability was again found in Zone 


IT along 18°E. The most striking feature is an elongated, 
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narrow zone of shallow thermocline depth, oriented NW to SE 
across most of the domain. This 1s believed to be an effect 
of the main branch of Atlantic Water flowing into the Ionian 
Sea. It corresponds reasonably well with the "tongues" of 
cool water found on the satellite IR images. In some areas, 
the main thermocline showed a step-like structure, yielding 
two different depths of maximum gradient. These areas were 
well-defined, but they seemed to be unrelated to any other 
particular thermal features. 

Surface ducts, due to a shallow isothermal layer, were 
present almost exclusively in Zone I. More interestingly, 
their temporal distribution was related to atmospheric forc- 
ing events. They appeared on June 9, persisted until June 14 
and reappeared for a short period of time on June 18 and 
June 20 (Fig. 36). Their average depth was 20m. 

A main sound channel was found all over the area, with 
an axis depth typically at 80 to 120m (Fig. 37). In subarea 
D, the distribution of the main sound channel was bimodal, 
with more than 50% of the axis depths shallower than 75m 
(see Table II). Subarea A had a deeper-than-average sound 
channel axis, about 130m. The thickness of the main sound 
channel varied from 800 to 1200m. Given the average bottom 
Semen Oc 3000m, the “depth excess" criterion [Urick, 1975, 
pl52] was met at almost every location in the Ionian Sea away 
from the shelf and the ridge, thus allowing the use of con- 
vergence zone detection. Ray tracing would indicate a 


convergence zone spacing of 34 km. 
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EaBLE If 


Frequency of occurrence of sound channel axis in each 
area. For the main channel, this frequency is indicated 
in percent for each layer, adding to a total of 1003 

for each area considered individually. The secondary 
channels are described by a rate of occurrence in each 
individual area, regardless of the depth at which it 

was found. Hence, it is an estimate of the probability 
of occurrence in that area. 


AREA 
A 8 eS 
No. of XBTs 14 50 26 20 49 
DEPTH ZONE MAIN SOUND CHANNELS 
35=75m 0 no 14 SS iO 
75-122m 23 54 46 44 80 
122-180m 77 20 Hes 0 10 
180-250m 0 Zi 27 0 0 


SECONDARY SOUND CHANNELS 


% occurrence 47 Ss S 0 5 


The common presence of temperature inversions below the 
thermocline is the key factor which explains the existence 
of secondary sound channels. Some of these small scale 
discontinuities were too tenuous to be accounted for, and the 
following limits were set, although somewhat arbitrarily, to 
filter the data: inversions were retained when their thick- 
ness was larger than 10m and the mean gradient larger than 
0-01 m/s m~'. In terms of sound propagation, the effect of 


such filtering was to reject sound channels usable at 
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frequencies higher than 5kHz (the low frequency cut-off for 

a 1Om-thick channel) and, at these high frequencies, diffrac- 
tive leakage losses higher than 13db/km (computed from 
Meourke, 1973}) would occur. Such restrictions would only 
apply to a limited number of operational situations. 

Secondary sound channels were almost nonexistent in 
Subareas D and E, which constitute Zone II, giving further 
evidence of the low variability there, as demonstrated previ- 
ously (Table II). The geographical distribution of the 
secondary features, overlayed on the main sound channel, 
provide further insight into the location of intense mixing 
mec, 37). 

The limitations of some equipments, in terms of accessible 
depth, might make the secondary sound channels appear attrac- 
tive as an alternative to the deeper main channel. Some of 
these secondary channels were continuous over extenSive 
distances (greater than 100 km), although many of them were 
only isolated (short-lived?) features. Several difficulties 
may arise when trying to insure the navigation of an acoustic 
array within such a channel, with an accuracy of about 5m in 
depth, in areas where intense vertical current shear might 
render the behavior of towed systems somewhat uncertain. 
Moreover, many factors could act on the thermal inversion to 
alter its characteristics: internal waves as well as inertial 


oscillations are the most common of these and have typical 
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periods of a few minutes to 20 hours (the local inertial 
period). The amplitude of vertical oscillations can be 
over 10m. Under such variable conditions, one might not be 
able to track the oscillations and keep an array in the 


selected sound channel. 
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iP eceMieeoON With THE BOTS ANALYSIS 


fe LOTS 

To allow some comparison between the present analysis and 
the EOTS (Expanded Ocean Thermal Structure) analysis, the 
following characteristics of EOTS are summarized. EOTS is a 
numerical scheme designed for the analysis of oceanographic 
data, mainly BT profiles and sea~surface temperatures. It 
is produced by FNOC, for use in the Fleet Oceanography 
Centers and forces at sea. It uses a climatological data 
base with both temperature and salinity as well as real-time 
XBTs. The Mediterranean Sea is covered by two analysis zones, 
Gach on a 63 by 63 grid, with a mean grid spacing of 40 km. 
The longitudinal boundary between the western Mediterranean 
analysis (MEDW) and the eastern analysis (MEDE) is approxi- 
mately 18°E. The analyses are run three times a week. Each 
input datum is weighted according to its nature and age, and 
blended with climatology. The most common outputs are 
horizontal charts of temperature at various levels, as well 
as the primary layer depth (PLD). Isotherms are contoured 
at 1°C intervals. 

The profiles used in the present study were among those 
used aS input to the EOTS, which may have used other inputs. 
The thermal analyses made in Zones I and II were compared 


with EOTS temperature analyses at the surface and 100, 200, 
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300, and 400m, obtained for June ll, 13, and 20 (Figs. 18-22 
(Zone I) and Figs. 30-32 (Zone II)). 

Not all figures display a synoptic scale feature, an 
indication that nothing was found or resolved by the EOTS 
analysis at that location. The deeper levels, 300 and 400m, 
in particular, were devoid of any features. When features 
are present in EOTS, they have an obviously coarser resolution, 
with a minimum radius of 40 km. This can be attributed to 
the numerical technique, which does not allow the resolution 
of features smaller than two grid lengths, i.e., 80 km in 
diameter. 

The present study, which spans an area from 14° to 19°E, 
required the use of both MEDW and MEDE analyses. It was 
observed that the two zones, where they overlapped, did not 
always show the same distribution of isotherms (Fig. 32). 
This was attributed to the fact that input profiles were 
selected geographically before each analysis, hence eliminat- 
ing the remote influence from data present on the other side 
of the boundary. 

Similarly, the primary layer depth (PLD) contours from 
EOTS were overlayed on the present thermocline analysis 
(Fig. 35). The comparison shows more consistency and better 
resolution than the temperature analyses. In particular, the 
limits of the elongated contours, which delineate a shallow 
PLD area, agree fairly well. However, the absolute values 


differ by 20m from one analysis to the other, the EOTS 
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showing shoaler values. This deviation is explained by the 
fact that the thermocline depth of this study is a contour 
of maximum temperature gradients, whereas EOTS computes the 
PLD as the depth of maximum gradient discontinuity, i.e., 


the top of the thermocline. 


eecolos INTERPOLATED SEA SURFACE TEMPERATURE 

Referring to Chapter VI (Figs. 18 and 19) and Chapter VII 
(Figs. 30 and 31), the EOTS interpolated SST was very helpful, 
and easier to contour than the XBT-derived SST. This demon- 
strates that at least this portion of the EOTS analysis 
schemes has excellent resolution and smoothing effects, 
allowing an accurate and consistent contouring. A quantita- 
tive comparison, which was not possible at depth due to the 
nature of the analysis grid, was made at the surface between 
EFOTS- and XBT-derived SSTs. The average difference was 
Semem> © (EOTS cooler than XBT), the standard deviation of 
this difference was 0.7°C, and the correlation of the two 
fields was 0.75. 

The author does not feel competent to analyze the numeri- 
cal processes through which the interpolated profiles are 
smoothed out and the information 1s lost when analysed on 
the 40 km grid scale. Whether it comes from grid averaging 
or from an excessive time-decay constant in the model is not 
relevant to this study. However, the quality of the inter- 
polation scheme is probably sufficient to justify an attempt 
to analyse the thermal field at a smaller grid scale, e.g., 


oor LO km. 
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X. COMPARISON WITH THE MILOC-68 INVESTIGATION 


MILOC-68 is the name of a NATO oceanographic investigation 
which took place in the Western Ionian Sea during the spring 
of 1968. It provided a number of STD casts as well as XBTs. 
The data reported by Miller [1972] were concentrated in a 
220x220 sq km area, during the 20-dav period from May 14 to 
June 2, 1968. This area represents the southern half of 
Zone II and the most northern part of Zone I (Fig. 5). The 
Similitude in both space and season with the present study 
was of great interest. In addition, the MILOC-68 report was 
received after the present analyses in Zone I and II were 
completed, which prevented the introduction of any subjective 


influence. 


Sweets MELOC-68 RESULTS AND COMPARISON 

The MILOC area spanned two very different subareas: the 
northern half was dynamically inactive while the southern 
half was dynamically active. The border between these two 
areas, lying along 36°N, was identified as a "high density 
Meemmal front” (Miller, 1972; p18]. The salinity and sound 
speed distributions exhibited the same baSic patterns (Fig. 
pea-d). 

A transect along 34°50'N showed the presence of a warm 
Seemaly (Fig. 39a). From a perpendicular transect along 


Pereeirigq. 39c), this anomaly was centered around 34°50'N, 
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18°E with a radius of 45 km, although some additional sampling 
south of 34°40'N would have been necessary to confirm the 
circular shape. 

AS pointed out by Miller, the depth of the 15°C isotherm 
is the most widely variable isotherm, a characteristic that 
has been found in the present study. The study by Saunders 
[1972] mentioned in Chapter V was made during MILOC-68, and 
is hence consistent and simultaneous with, if not identical 
to Miller's analysis. The freshest water was in the south- 
western part of the area (Fig. 38e) and correlates remarkably 
well with the tongue of surface water, 1.e., Atlantic Water, 
flowing eastward. To the north, high salinity LIW "encroach- 
ed upon the shelf at the entrance to the Sicily passage 
between Malta and Sicily" [Miller, 1972] as an identifiable 
water mass. Both temperature and salinity fields can be 
related to the present study, and the characteristic 21°C 
Surface isotherm (Figs. 30 and 31) can be taken as a limit 
between the stationary warm area (north of 36°10'N), 
presumably with a higher surface salinity, and the more 
perturbed area. To the south, the mixing area showed the 
following common features: 

the cool tongue seen on satellite pictures (see Figs. 
14 and 38a-e) 
mm@emcool tongue (20°C isotherm) found in Zone II (36°N 


17°E) (see Figs. 30 and 38a-e) 
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the warm eddy found in Zone I (34°30'N 18°30'E) and 
Titec esttmenat tO ehemMmnOe-o38 “eddy” (34°50°N 18°E) 


(see Figs. 24 and 39). 


ferib INTEREST OF THE MILOC-68 STUDY 
From the comparison made previously, it can be stated 
maa t : 
- there is a considerable consistency between the two 
studies, 
the present study was more sketchy, due to the random 
sampling, the non-synopticity in Zone II, and the lack 
of salinity data. 
A data base such as MILOC-68 is very valuable in many respects. 
First, it can provide the analyst with an initial impression 
on which he can build, by successive approximations, a con- 
temporary description of the area, provided the seasonal and 
Spatial domains and oceanic and atmospheric conditions are 
Similar. Second, it provides a high quality case study, more 
carefully checked than is usually done operationally with 
mee Finally, it has salinity profiles which, although not 
essential to sound speed computations, are important to any 


investigation of dynamical processes. 


eae PLICATIONS 
Given the similitudes described above, the MILOC-68 data 
were used for some additional analyses, assuming they applied 


reasonably well to the present study. The first application 
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attempted to relate the T-S characteristics of 13 CTD profiles 
to their surface and subsurface structure. The second appli- 
Cation was a computation of baroclinic modes and the Rossby 
radius of deformation. 

1. Correlation between surface and subsurface structure 

The 14 profiles present in the 34°50'N transect (Fig. 
39a) were analysed on a TS diagram (Fig. 40). The profiles 
were numbered from 1 to 14 from west to east, profiles 9 to 
12 being the warm anomaly. (Profile 14 was withheld from the 
analysis due to a suspicious anomaly in salinity at the 70m 
level.) 

At the surface, there is a dense grouping of profiles 
mmeeemo (Group 1) around 20°C, 37.6%. Similarly, profiles 
Meee !2 (Group 2) are concentrated around 21°C, 37.9%. 
Profile 13 1s singular, at the same temperature as Group l 
but with a somewhate lesser salinity. 

Between 50 and 100m, almost the same grouping is 
present, but profiles 1 and 3 are in Group 2, which is still 
more saline and warmer than Group 1. Around 100m, all pro- 
files merge towards the local LIW water type at 14.6°C, 
38.83%. 

Succ GEOUpIng Of profiles 1s not random; in fact, 
it can be associated closely with the structure found in the 
transect and provide some knowledge about the composition of 
the upper part of the warm anomaly. The profiles in Group l, 


in the 50 to 100m layer, are indicative of mixing between 
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surficial Atlantic Water and local LIW. By contrast, Group 
2 is representative of mixing of local LIW with a warmer, 
more saline surficial water type. This is characteristic of 
surface water found to the east. The upper layer of the 
warm anomaly could have been formed in a mixing event to the 
east, a hypothesis which suggests a westward drift of this 
warm anomaly as an isolated dynamical entity. The layer 
considered here is below the thermocline, and a dynamical 
process (overturning, mixing, etc.) is more likely to have 
produced it than penetration of heat across the thermocline. 

Furthermore, profiles 9 and 12, located on the edges 
of the warm anomaly, were closer to Group 1 than were profiles 
10 and ll. These two profiles, being located in the core 
of the anomaly, could possibly have suffered less mixing 
than profiles 9 and 12. Profile 13 was even less saline 
than Group 1, which indicates the presence of a surficial 
lens of fresher water, for which no physical explanation is 
advanced. 

To explain the behavior of profiles 1 and 3, which 
demonstrated characteristics similar to those of profiles 9 
to 12 at the 50 to 100m level, one can refer to the tempera- 
ture and salinity transects which indicate the edge of a 
trough-like perturbation. It is accordingly hypothesized 
that they might represent the same water mass composition 
eeeprotiles 8 to 12, but did not show a tendency toward a 
"Group 2" composition at the surface because of some 


B@inricial transformation. 
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2. Internal (baroclinic) wave computation 
"Internal waves are fluid motions which are under 


the influence of variations in the ocean density (buoyancy) 


field and the earth's rotation; i.e., both these influences 
provide restoring forces" [Mooers 1975, pl067]. The baro- 
clinic low-frequency wave propagation in a rotating fluid 


is governed by the following equation: 


N*(z) W + f° W = 0 
xo Ze 


with o*<<f* where W(x,z) represents the vertical velocity, 
N(z) the Brunt-Vaisala frequency, o the wave frequency, and 
f the Coriolis parameter. Assuming that W(x,z) can be 
modeled as: 

niece 2) = 4(zZ) exp (kx), 
the previous equation becomes: 


eo Keates — 0 


This eigenvalue equation was solved for k with the following 


kinematic boundary conditions (BC): 


mij y= 2(=D) = 0 





i.e. the vertical velocity vanishes at the surface and 
bottom. A Brunt-Vaisala frequency profile was computed from 
profile 11 (MILOC-68), the center of the warm eddy. Once 
the problem was solved for the first five modes, the eigen- 
functions Z(z) were evaluated (Fig. 41) and the baroclinic 
Rossby radius of deformation computed as R = 1/k (Table III). 
Values of Rare commonly interpreted as typical scaling 
lengths for dynamical processes. The eigenfunctions also 
give some information about the vertical dependence of 
velocity fields, and show in particular the depths where 
maximum shear and horizontal velocity reversal occur. 
Details of the computational algorithm are given in 


Appendix B. 


TABLE III 


Values of the baroclinic Rossby radius of deformation 
(modes 1 to 5). 


Mode Number IL 2 3 4 5 


Rossby radius (km) 26 2 9 a 5 


Assuming that most of the energy of an eddy is 
associated with the barotropic and the first baroclinic mode, 
23 km is, hence, assumed to be a reasonable first-order 
scaling length for the perturbation field. This value can 


be related to the e-folding distance of the correlation 
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Pmeigiettons touma im Zone - and II, 1.e., 30 and 45 kn, 
respectively. However, a visual comparison of the second- 
mode vertical structure with the isotherm perturbation field 
(Fig. 39) shows that the main 15°C isotherm perturbation 
occurs at the same depth as the maximum of the second mode 
vertical speed. Hence, the second mode may also be important 
to describe the spatial structure, with a Rossby radius of 
deformation (second mode) of 12 km. Given the inherent 
approximations, including the fact that the Rossby radius of 
deformation was computed from a totally different data set 

at the center of the most notable anomaly, and with no data 
below 400m, the agreement is fair. The question is raised 

as to whether the local computation of the Rossby radius of 
deformation” could be used as a valuable operational index of 
moerizontal variability. 

The development above was limited to the f-plane, 
mew GO an ocean with no latitudinal change of the Coriolis 
parameter. On the 8-plane, where f is allowed to vary with 
latitude (8 = 3f/3y), the theory indicates that a class of 
long waves, called planetary or Rossby waves, is allowed to 
propagate in the (baroclinic) ocean. The wavelength (L) of 
such waves is of the order of 27k, the period (T) is greater 


than the inertial period: 
[R7*+ 4] 
T= (27)? ———__ 
8L 'cosé 


ore actually treating a baroclinic wave, the avprovriate 


scale is the horizontal wavelength, which is of the order 
of 2TR, or about 72 km in the case of a second mode wave. 
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where 98 is the direction of propagation relative to east and 
the phase speed has a typical order of magnitude of a few 
cm/s for an internal Rossby radius of deformation of about 


10 km, which always has a westward component. 
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XI. OCEANOGRAPHIC SUMMARY 


The various analyses reported in the previous chapters 
led to a description of some oceanic characteristics of the 
Ionian Sea. None of these characteristics seemed to differ 
from the commonly accepted picture of the region, but some 
of them bring a clearer focus to specific points, and 
suggest new questions. 

The series of satellite IR images provided some quantita- 
tive estimates of the temporal variability and the geographical 
extent of the surface advection patterns, as the surficial 
water flows from the narrow Strait of Sicily into the more 
open Ionian Sea. Questions arise regarding the factors that 
influence and determine the characteristic double-tongued 
flow: where is the generation area and what are the trig- 
gering factors for this splitting? Is there a vortex pair 
(a cyclonic and an anticyclonic eddy) formed? What processes 
Bevermm the anticyclonic curling of the northern tongue? What 
is the influence of bottom topography; e.g., the continental 
Shelf and the ridge? Once this flow pattern is confirmed 
and its seasonal variability evaluated, a laboratory simula- 
tion could be undertaken to assess some of these factors, in 
a way Similar to that which has been done recently for the 
Petatt Of Gibraltar {Whitehead and Miller, 1979]. A dynamical 


(numerical) approach could also be attempted to verify, for 
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each of the water tongues, whether the flow pattern could 
be affected by the release of absolute potential vorticity 
(baroclinic instability) as it reaches the deeper ocean 
[Mysak, 1977]. 

The satellite IR imagery correlated well with two sub- 
surface features. First, the distribution of the thermocline 
depth had a good visual correlation with the advection pat- 
terns of surficial Atlantic Water. A relatively shoal 
thermocline depth can be expected in these tongues of cooler 
water. Second, the presence of a warm anomaly was revealed 
by the satellite IR imagery, as a transient atmospheric 
boundary layer effect, with the occurrence of low level 
Stratus. Boundary layer theory indicates that such events 
would be most likely in the early morning and in the presence 
of higher winds. The observation of such a feature, if 
correctly interpreted as an atmospheric "second-order" effect, 
can be a first suggestion of the presence of such a feature, 
and may allow its (discontinuous) tracking. 

The XBT analysis revealed - or confirmed - some contrast 
between the northern zone (Zone II) and the southern zone 
(Zone I). On one hand, Zone II is a relatively quiescent 
area in which the superposition of water masses has no 
dramatic consequences, and little apparent (Synoptic scale) 
Memeezoncal variability; on the other hand, Zone I is directly 
exposed to the entering surface flow from the Strait of Sicily, 


mmrcnh seems to create geophysical turbulence and baroclinic 
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instability, as revealed by the occurrence of temperature 
inversions and layering. One particular feature was a warm 
core eddy, about 60 km in diameter, characterized by a strong 
(200m) deepening of the 15°C isotherm. A warm core eddy was 
also found in another study (MILOC-68); hence, it suggests 
that eddies could characterize the synoptic scale turbulence 
in that region, although with possible seasonal variation. 
Questions arise regarding the origin of this warm anomaly, 
which has a different water mass composition from the ambient 
water. It was suggested that this anomaly could have been 
formed by mixing to the east of its location. This hypothesis 
would rewuire further investigation to explain the westward 
migration of this isolated feature, perhaps as a Rossby wave, 
or eddy, and its associatin with depressed isotherms at the 
200m level. Such an investigation could be possible, at first 
approximation, using a two-layer’ model of the ocean and 
disregarding any bottom topographic effect [Leblond and Mysak, 
1978] or using a more sophisticated approach [McWilliams and 
Flierl, 1979]. Scales for wavelength and phase speed could 
be estimated for comparison with observations as in 


maeurgut [1976]. 


Te the second mode structure noted earlier is predominant, 


a three-layer model would be needed. 


78 





Finally, the statistical study of spatial and temporal 
scales of variability showed a clear anisotropy in Zone I. 
How typical these correlation functions are, and how perma- 
nent they are over an annual cycle could not be assessed, 
and would require a larger data base. It is suggested that 
such two-dimensional spatial-lag correlation functions would 
be useful to better define the spatial scales required for 
sampling and analysis. Clearly, additional analytical 
techniques are needed to allow for a systematic separation 
of the different processes that distort the synopticity of 


the data, which hence, decrease the quality of the analyses. 
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XII. OPERATIONAL CONSIDERATIONS 


me LWO OBJECTIVES FOR OPERATIONAL ANALYSIS 

The basic methodology used in the present study can be 
described as an attempt to analyse and interpret a temperature 
data set in the light of the regional oceanology. It is the 
belief of the author that such an approach should be used as 
a routine technique to provide navies with environmental 
ireaucts. 

It is a usual procedure for a group of naval ships to 
organize a bathy(thermograph) watch by assigning a guard ship 
with the task of sampling at regular intervals, usually 4 or 
6 hrs, and broadcast each profile to other units, as well as 
to the environmental facilities ashore. Units then process 
this single profile to obtain the variables they need: 
predicted range, optimum depth, etc., for the appropriate 
sensors. In some tactical situations, the spacing between 
units might be re-evaluated as well. 

But the sampling is seldom organized according to the 
@eean itself, its variability, the atmospheric forcing 
processes or the location of previous data samples in the 
region. Moreover, airplanes and helicopters have their own 
procedures, which thev most often follow independently of 
M@ensurcace units. In fact, the error in such an approach 


is to reduce a four-dimensional problem (time, space) to a 
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two-dimensional one (time, depth). This procedure has two 
major deficiencies: 

First, it makes sampling a random operation, with a loss 
of information when areas of strong variability are crossed 
and not sampled, and an increased expenditure when XBTs are 
redundant in unstructured areas. 

Second, it encourages ASW forces to regard the ocean as 
a horizontally homogeneous water mass, with a temperature 
structure varying incrementally every 4 or 6 hours. It most 
certainly causes a lack of accuracy in sonar range predictions. 
But worse, it precludes the ability to use such inhomogeneities 
in the ocean for tactical purposes: hiding a unit, optimizing 
a silent transit, enhancing detection capabilities, etc. 

Ashore, prediction facilities do not have the appropriate 
tools to conduct a more thorough analysis of the ocean 
structure, at a scale compatible with operational requirements. 
For example, at FNOC computational techniques are currently 
used which blend a variety of input variables through weight- 
ing of these variables [Holl et al., 1979] over extensive areas. 
This approach might be of some interest as a method of producing 
an input to global atmospheric models. However, it does not 
seem to be, as yet, well adapted to meet naval operational 
ocean prediction requirements, which involve typical scales 
of 5 to 100 km, suggesting a need to resolve, at least, the 


20 to 60 km scale, so commonly found in the ocean. 
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Finally, reasonable objectives for ASW forces in this 


attempt to better understand and exploit the ocean would be: 


i. 


to routinely include the horizontal variability, through 
a more appropriate sampling procedure, and 

to request synoptic scale analyses from shore prediction 
facilities, which implies a regional approach to 


oceanography. 


These two points are addressed in the following sections, 


with the idea that the next step in this development would 


be the implementation of short term, on-board prediction 


tools. 


5. 


IMPROVING DATA COLLECTION 


Four complementary systems 


A variety of systems are available currently for 


monitoring the oceanic structure in a region: 


XBT profiles are well known for their ability to 


describe the temperature and, hence, sound speed 
structure, in the upper layers. This study has shown 
a potential capability to evaluate the horizontal 


fields. 


Satellite remote-sensing is available to provide a 


variety of variables, among which SST is the most 
common. The main interest in such data is to provide 
the analyst with a perception of synoptic scale and 


mesoscale features, and their approximate location. 
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° Drifting buoys, even a Single one, could bring some 
dynamical clues in the evaluation of time variability 
(heating and turbulence of the upper layer), as well 
as advection, by continuous tracking. 

- Continuous SST (and salinity) measurements from surface 
Ships, if accurate, could clearly indicate anomalies 
such as fronts or eddies, and aid in the interpretation and 
"Calibration" of satellite IR images. Furthermore, 
they could better orient XBT sampling in real time and 
encourage a better perception of horizontal variability. 

2. Recommendations for XBT sampling 

It 1S not known what type of XBT sampling procedures 
were implemented during the SHAREM 38 exercise from which the 
present data set was acquired. The limits of the exercise 
zone were also not known. Hence, no appreciation could be 
gained as to how well the actual XBT sampling covered the 
zone. In fact, the sampled area might have been somewhat 
smaller than the operating area, a situation which commonly 
occurs due to obvious priorities and limitations in 
Meaility. 

Fifty-five XBTs constitute the coverage in Zone I, 
an area of 30,000 sq km. Some pairs of XBT profiles were 
positioned quite close, in both time and space, and could be 
termed as a redundancy in the data set. This redundancy 
could be prevented and yield a better coverage within the 


mmres Of a planned sampling effort. 
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The correlation diagrams presented in Chapter VI 
suggest an order of magnitude for spacing of XBTs, Since it 
evaluates the amount of common information present in two 
profiles. In particular, 1t suggests very strongly that 
Sampling should not be isotropic. Since the Ionian Sea is 
clearly much more variable in the E-W direction, sampling 
Should be made at closer intervals in this direction. The 
+0.5 correlation contour could be used as a Satisfactory 
level of cross-correlation, suggesting that 20 km in the E-W 
direction might be adequate (Figs. 27 to 29). Spacing in 
the N-S direction depends on the layer of interest: adequate 
Spacing would be from 30 to 60 km. 

This type of approach, where correlation functions 
in a layer are used to estimate adequate sampling intervals, 
Should be consistent with requirements for resolving features 
in the sound speed field. This 1S a factor which depends on 
sensor characteristics, mainly frequency and depth. Assuming 
that the depth of minimum sound speed, which is between 100 
and 200m in this region, might be of interest, a sampling 
distance of 60 km would be adequate in the N-S direction 
(Fig. 29). Under such a procedure, Zone I could have been 
adequately covered with 26 XBTs. In other words, the same 
XBT expenditure would have permitted sampling the entire 
area twice within the same five-day period. The term 
"adequate coverage" implies here an attempt to resolve 


Synoptic scale features, in order to track them and take 
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them into account in naval operational planning. This does 
not suggest that XBT sampling at short (a few hours) intervals 
is irrelevant, especially if IR imagery is used to guide the 
pattern of XBT deployments, but it emphazises the interest 

of well planned spatial sampling and analysis, by contrast 
with a temporally regular/spatially random procedure. [In 

the upper layer of interest to hull-mounted sonars, more 
attention must be devoted to the transient physical processes 
and closer sampling will hence be required, until the driving 
processes are better understood and forecasted. 

Clearly, the analysis grid, just like the sampling 
grid, should be determined carefully with respect to the 
thermal structures that one tries to resolve. The present 
study used the sampling grid as an analysis grid, since the 
data were not overwhelming, and the thermal fields were evalu- 
ated and smoothed manually. It is nevertheless interesting to 
proceed faster, and more objectively, i.e., consistently, with 
computer assistance, and objective analysis techniques (using 
weighting factors based on the correlation functions) are cur- 


rently developed for such purposes [Bretherton et al, 1976]. 


C. AN APPROACH TO REGIONAL OCEANOGRAPHY 

The suggestions developed in the following section apply 
to a geographical region where continuous monitoring of the 
oceanic variables is required, whether it is for harbor 
defense, the control of vital shipping lanes or the surveil- 


lance of submarine traffic. It could be adapted to the 
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transit of a task force, which usually covers a long track 
through a variety of oceanic regimes. 

First, areas of operational interest should be delineated, 
and their boundaries should be kept consistent with the 
oceanic regime and local processes. 

Second, a compilation of experimental studies should be 
made for each area, to provide the oceanographer with the 
initial material, a step which is instrumental in the region- 
al approach. Such a study should be physical as well as 
statistical and describe processes such as formation of deep 
water, circulation, mixing of water masses, typical variables 
or products representative of these processes, etc. Atmos- 
pheric forcing variables should not be overlooked, and in 
fact need to be closely associated, especially with respect 
to mixed layer processes, Ekman transports, and formation of 
deep water. Such a compilation of local information can be 
compared, at a different scale, with forecasters' handbooks 
that are developed by meteorologists to provide coverage of 
local conditions, particularly near airports [Brady and 
Nestor, 1980]. 

Third, the sampling should be organized and adapted to 
mae rbegion, 1.@€., to both spatial and temporal scales of 
dominant processes. This point has been more fully developed 
in the previous section. 

Fourth, a set of standard tools should be made available 


to help the analyst to process incoming data and relate them 
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Pomcne climatological picture. The evaluation of how far 

the computer software must be involved in the data processing 
is precisely one of the steping stones of this method. These 
programs should not be used beyond the point where they start 
removing the essential information from the input data, 
leaving only the climatological picture. In fact, some 
computer-assisted manual evaluation is very likely to remain 
Meerul for a long time, which applies to quality control, 
field contouring and any other non-mechanical technique, a 
fact which is acknowledged in their own field by many 
experienced meteorologists. 

Lastly, more specialized analytical tools should be 
developed as research progresses, and could be expected in 
the near future in such fields as data integration, dynamics 
of frontal regions, Rossby wave propagation, or mixed layer 
evolution. There would be the main tools necessary for 


operational analysis and forecasting. 


D. THE MEDITERRANEAN SEA, NATURAL LABORATORY 
The Mediterranean Sea offers an interesting potential as 
a natural laboratory for coupled environmental/tactical 
experiments: 
1. Being surrounded by continental masses, its synoptic 
weather situation is well monitored and does not suffer 


much from data gaps. 





Peelosegrapareal peculiarities like straits, sills or 
baSins are abundant, and create a variety of regimes 
in both the regional and dynamical sense. 

3. The conditions for the use of satellite IR remote- 
sensing for the determination of thermal properties 
and other surface variables are exceptional, due to 
low cloud coverage throughout most of the year. 

4. The "moderate" size of its basins, and their clear 
physical separation could allow for synoptic-scale, 
three-dimensional, numerical modeling which is within 
the reach of current computer capabilities once the 
dominant physical processes are properly determined. 

Given these natural qualities, the Mediterranean Sea 
appears as an ideal region where new oceanographic techniques 
and tactics could be tested and improved, prior to their 
application, mutatis mutandis, to other regions of ocean, 


where the conditions are less exceptional. 
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APPENDIX A 


oleolueat et orics 


ee) DEFINITION 


The following notation is defined: 


n the total number of synoptic XBTsS in a zone 
1,j the indicies defining particular profiles 


k the index for depth for a given profile 


Hence TZ; ke represents the temperature at level k, l.e., 
ti 


at depth 5(k-1) meters from profile at position 1. 


A.2 METHOD 

Each profile was digitized by 5m increments from surface 
(k = 1) to 400m depth (k = 81). A mean profile was then 
computed by averaging the temperature at each digitized depth 
over the whole set of n XBTs: 
Siven a Dalr Of XBT profiles, a cross-correlation 


coefficient was computed for each previously defined 


layer: 
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in the 0 to 60m layer; and similarly in the 60 to 150m layer 
(with summation from k = 13 to 31) and the 150 to 400m layer 
(with summation from k = 31 to 8l). 

Given a set of n profiles, the number of cross-correlation 
coefficients generated in the process is n(n-1)/2 . The 
geographical location of each profile in a pair yields a ARG ss 
eer... 


Bs 


these two profiles. Whenever AX, 5 is found to be negative, 


both signs of Ax, 5 and AY 5 must be changed. This is essen- 


defining the horizontal E-W and N-S distance between 


tially the same as limiting the azimuth between XBTs to the 
O0-180° range without loss of generality. 


AY...) were related to the 


mae set of values (p.. iy’ i3 


eg 
time interval in the following fashion: the distance AD; 
between the two profiles was computed and associated with 
the difference in days between the two profile dates. Dis- 
tances were quantified and averaged by 10 km increments, 


whereas time intervals were expressed in an integral number 


of days, from zero to four (Zone I) or seven (Zone II). The 


oe 





cross-correlation coefficients were plotted in this time- 
Space coordinate system. 
AX A 


Next the sets of values ( ) were quantified 


oe Oa i CT 
by 10 km intervals, and values of the cross-correlation coef- 
ficient falling in a given 10 by 10 km Square were averaged, 
yielding a single coefficient. 

Finally, these averages were plotted on a 10 km square 
grid and contoured, for each of three layers. The number of 
cross-correlation values in each square was also plotted to 
give a measure of the density of information. Clearly, these 
plots are symmetric with respect to the origin, since it does 
not make any difference whether profile 1 is considered to be 


memene center, and profile j in the AX, AY grid, or the 


Seenverse. 





APPENDIX B 


NUMERICAL CALCULATION OF BAROCLINIC INTERNAL WAVES 


B.l1 EIGENVALUE PROBLEM 


Solving equation (1), with its associated boundary 


conditions (BC), is an eigenvalue problem: 
meet N  (z) k* £-°2(z) = 0 ip) 
Z{0) = Z{-D) = 0 


Pwesolucion for k and Z can be found for an infinite number 

of modes. N(z) is a geophysical variable, the Brunt-Vaisala 
frequency, which represents the natural frequency of oscilla- 
tion of a water particle under the influence of the buoyancy 


force: 


30 a0 2 2 
N*(z) = -g/= (32 = -g/- (32) = Ca 7e e 
e Oz adiab p 3Z ims ae ee 


where o0(z) is the density of water, g the gravitational 


acceleration and c(z) the sound speed at depth z. 


B.2 NUMERICAL METHOD 
For the Brunt-Vaisala frequency profile N(z) encountered, 


an analytical solution to equation (1) 1s not possible. A 


‘Oo 
Or 





numerical technique obtained from Bell [1971] was used as 
follows: 

The Brunt-Vaisala profile is digitized from the surface to 
the bottom, yielding n values Noe Noe ---N, at n levels 


Zs ee 2262). These levels need not be at evenly spaced 


intervals. An approximate quadratic form for Z(z) 1s assumed 


for each depth interval spanning two increments of depth. 


Z =a + a ; +a Z 

es | o1 Ae) a) 2a Ale ot 
7, =a ta. 2a + A Zo 

ale GL Neale 1 9 i. E 

Ae =e ee a. 3, oh Alok os 
et] jal 11 let a1 iLse) 


The second derivative Z''(z) 1s equal to 2aiae and the system 


of equations above can be solved for a as a linear function 


mat 


Qs : as = O(Z5_, air eee 


Hence, equation (1) becomes: 


meee NCE *Z. = O 
a a 


Or 
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This can be solved for Z. 
Lt] 


ae aio aD 
: = — — : = : =_ : . Ns : 
are nea t [=k Nit ae 2 ee 43] ay - lida > ree 
(2) 
where A. = 25,4, 72; 
2 
and 1 Ci ee 
det eae ec 
ak 
4 Si 1+; 


eee lhe LTERATION SCHEME: A SHOOTING METHOD 
The theory of eigenvalue (Sturm-Liouville) problems 
asserts that the wave number which yields solutions satisfying 


the BCs must lie in the interval: 
Ome Kae eee ee ole 
; | 


Also, the first mode eigenfunction has one and only one 
extremum between the boundaries, hence no zero-crossing, 
the second mode has two extrema and one zero-crossing and 
Zemeorth. 

The shooting method for the first mode is designed as 
Follows: An initial value for k is guessed; given a value 


eo ay =O Sur tace 3e) and Z = 1] (normalized value), the 


3 


Z.'s are computed from equation (2) by iteration to the nth 
iMevel, 1.€. the bottom. If the axis Z = 0 is crossed before 
reaching the nth level, a new lower value for K 1S ae ered . 
Otherwise, the value of Z 1s examined, as the bottom BC. A 
larger value is tried for k until Z, gets reasonably close 
to zero. When the iteration is stopped, k is the numerical 
solution for the first mode wave number, and the set of Z.'Ss 
(i = l,n) is the eigenfunction. 

The next mode iS initiated with a guess value K = 2k - 
The iterative process is conducted to allow for the proper 
number of axis crossings, and to obtain a value Z = 0; and so 
forth for the other modes. The eigenfunctions obtained are 
also representative of the density eigenfunctions, and the 
general theory of baroclinic waves show that the horizontal 
velocity and the pressure eigenfunctions are proportional to 
the first vertical derivative of the eigenfunction Z [LeBlond 
and Mysak, 1978]. These vertical profiles were computed from 


each eigenfunction Z by a simple centered differentiation 


scheme. 
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